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Abstract 

■g : A new type of cloak is discussed: one that gives all cloaked objects the appearance of a flat 

I conducting sheet. It has the advantage that none of the parameters of the cloak is singular and can 

P3 ■ in fact be made isotropic. It makes broadband cloaking in the optical frequencies one step closer 

I— J 

\0 ■ PACS numbers: 42.79.-e, 02.40.-k, 41.20.-q 
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Transformation optics l|, |2|, y, |J] can be used to design a cloak of invisibility. In essence 
the cloak makes its contents appear to be very small and hence invisible. However there 
are three distinct topological possibilities: the cloaked object can be crushed to a point, to 
a line, or to a sheet. In the process of crushing the object becomes infinitely conducting 
but this does not present a problem for the first case because scattering from even highly 
conducting small objects vanishes with the size of object. The same is true of a very thin 
object, i.e. a wire, because very thin wires have very large inductance and are therefore 
invisible to radiation. This is not true of the third possibility. Obviously a conducting sheet 
is highly visible unless of course it sits on another conducting sheet (the 'carpet' of the title). 
Although this third possibility has limited cloaking potential it does have the considerable 
advantage that the parameters of the cloak need not be singular and, as we shall show, can 
be isotropic. 

Metamaterials are usually used as the building blocks for transformation optics as they 
can attain a wide range of permittivity and permeability. Several works investigated the 
cylindrical cloak in full wave analysis and the first experiment in verifying the cloak is done 
m microwave. U U m m There are other approaches targeting the subwavelength regime 
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[91. II OH and the geometrical optics regime 11| for cloaking as well. Here, we concentrate on 
the transformation optics approach which is driven by coordinate transformation and has a 
large variety of applications. Q 

There is still a general quest to push the working frequency to visible regime for more 
applications. However, scaling down the magnetic resona ting structures for the metamate- 



rials to optical frequencies causes severe absorption. llSl . Il6l . llTI . Il8| One way to avoid 
it is to use the reduced material parameters in the cloak so that only anisotropic metama- 
terials with electric resonating elements are needed, [l^ However, as long as we are using 
metamaterials with resonating structures, materials absorption still cannot be neglected. 

In this work, instead of the complete cloak, we will consider a cloak to mimic a flat ground 
plane. We will show that it does not require singular values for the material parameters, i.e. 
the range for the permittivity and the permeability is much smaller than in the case of a 
complete cloak. Moreover, by choosing a suitable coordinate transform, the anisotropy of the 
cloak can be minimized to a small value also. As a result, we can avoid using metamaterials 
and just use isotropic dielectrics to construct the cloak. It greatly reduces absorption and 
also make broadband cloaking one step closer. 
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Here, we consider only the 2D wave problem for simplicity, i.e. all the fields are invariant 
in the z-direction. Furthermore, the E-polarization is assumed while the formulation can 
be easily adopted for the other polarization. A ground plane here means a highly refiecting 
surface made of metal. It is regarded as a perfect conductor. Suppose an object lies on 
it, we want to design a cloak covering on the object so that the observer perceives the 
system as a fiat ground plane again without introducing any additional scattering from the 
object. The object is now concealed between the cloak and the original ground plane as 
shown on the left hand side of Fig. [1] We further assume that the cloak (shown in cyan 
color) has a rectangular shape of a width of w and a height of h except that the bottom 
(inner) cloak boundary is curved upwards to leave enough space to conceal the object. 
The whole configuration is termed as the physical system with coordinate (x, y) or literally 
written as (x^, x^) in indexed notation. The virtual system is the configuration the observer 
perceives. It is shown on the right hand side of Fig. [H Its coordinate is labeled by (^, rj) or 
literally written as (^^,^^) in indexed notation. In general, we are considering a coordinate 
transform which maps (with sense preserving) a rectangular region (0 < < w,0 < r] < h) 
in the virtual system to an arbitrary region (the cloak in this case) in the physical system. 
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FIG. 1: (color online). The general mapping between the virtual system and the physical system. 
The regions in cyan color are transformed into each other. The shaded region at the bottom 
of both domain represents the ground plane (a perfect conductor). The observer perceives the 
physical system (x-system) as the virtual system with a flat ground plane. 

To ease our discussion, we introduce the Jacobian matrix A by 
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A* 



w 



(1) 



and the covariant metric g by 



9i'j' 



el-g or g = A^A, 



(2) 



where ^i, ^2 are the covariant basis vectors of the virtual coordinates appearing in the phys- 
ical system. In the application of cloaking, the observer perceives the physical system as 
the virtual system of an isotropic homogeneous medium of permittivity e^gf and unit per- 
meability. The corresponding permittivity and permeability in the physical system induced 
by the coordinate transformation are given by 



e = Eref/ Vdetg, 



y/det g 



AA^. 



(3) 



We can write and /i/, be the principal values of the matrix representation of the perme- 
ability tensor in the physical domain and the corresponding refractive indices be ut = y/f^i^ 
and riL = \/\i'T^ for the two (local) plane waves traveling along the two principal axes. To 
indicate the extent of anisotropy of the physical medium, the anisotropy factor a (a function 
of position) is defined by 



a = max {nx/nL, ni/nT) 
By using Eq. ([3]) and Eq. ([2]), it can be proved that 



(4) 
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Tr (g) 
V det g ' 



(5) 



with 



(6) 



On the other hand, we can define an averaged refractive index n relative to the reference 
medium by 



(7) 
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so that 



Sref yjaetg 

Instead of using e and /i*-^ to describe the physical medium, we now use a and n (related 
to the ratio and the product of ut and ul) which have geometrical meanings in terms of the 
metric. If we imagine there is a very fine rectangular grid in the virtual domain with every 
tiny cell being a square, the mapping (or the coordinate transform) transforms this grid to 
another grid in the physical domain. Every such tiny square {5 x 5) in the virtual domain 
is transformed to a parallelogram with two sides ^i5and .^5. A smaller anisotropy means a 
smaller value of Tr (g) / ^/detg while a smaller area of the transformed cell {\/det gS^) means 
a larger refractive index n. 

In cloaking, compression of space in the physical domain essentially makes the cloak 
anisotropic. However, our heuristic approach is to minimize the anisotropy induced in the 
physical medium by choosing a suitable coordinate transform. If the anisotropy is small 
enough, we can simply drop this part (by assigning a = 1) and only keep the refractive 
index n. In other words, the physical medium becomes just a dielectric profile described by 
Eq. ([8]) with unit magnetic permeability. 



In t 
tional 



lis work, such an optimal map is generated by minimizing the Modified-Liao func- 
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I ru> rh Tr{gf 
hw Jo Jq det g ' 
upon slipping boundary condition. Slipping boundary condition means that each of the four 
bounding edges of the virtual domain must be mapped to the four specified boundaries in 
the physical domain, only up to a sliding freedom. The minimal of this functional occurs at 
the quasiconformal map 2l| . Without going further for technical proof, we state the result 
that the quasiconformal map actually minimizes not only the average but also the maximum 



value of Tr (g) / ^/det g in the physical domain, i.e. anisotropy minimized. If we imagine 
there is a very fine rectangular grid in the virtual domain with every tiny cell being a square, 
the transformed grid in the physical domain for the quasiconformal map has the property 
that every transformed cell is a rectangle of a constant aspect ratio M : m where M is 
called the conformal module of the physical domain (a geometrical property of the physical 
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domain once the four boundaries are specified) and m = w/h is the conformal module of 
the virtual domain, i.e. 
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m 



a/ det J 



By substituting Eq. (ITO!) into Eq. ([5]), we have 



(10) 



Tr (g) M m [Mm 

, = h — , or a = max — , — 

^/detg m M \m M 



independent of position 
1.5 




(11) 



FIG. 2: (color online). The transformed grid in physical system with inner cloak boundary specified 
by Eq. (jl2p for the (a) transfinite grid and the (b) quasiconformal grid. The color maps show 
profile n?. 

As a first example, let us map the area of a rectangle bounded by < < 4, < < 1.5 
in the virtual domain to the same rectangle in the physical domain but with the bottom 
boundary specified by 



^bottom 



XI 



0.2cos(7rx/2)2 1 < x < 3, 
otherwise. 



(12) 
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The transformed area in the physical domain is exactly the region of the cloak. Fig. [2]^a) 
shows the transformed grid in the physical system if a simple transfinite interpolation is 
used to map a regular 40 x 15 grid in the virtual domain. In this case, the grid is just a 
linear compression in the y-direction. The corresponding range of the anisotropy factor a, 
obtained from Eq. ranges from 1 to 1.385 while n^, obtained from Eq. ([H]), ranges from 
1.0 to 1.153. On the other hand, for the generated grid using quasiconformal map, shown 
in Fig. [2t^b), the grid lines are orthogonal to each other. The aspect ratio of each cell or the 
anisotropy factor a becomes a constant of 1.042 while ranges from 0.68 to 1.96. As the 
area of each cell is proportional to the square reciprocal of n, the minimum of n occurs on 
the inner cloak boundary at around x = 1 where the cell area is largest while the maximum 
occurs also on the inner cloak boundary at a; = 2 where the cell area is smallest. Therefore, 
we sacrifice n to a larger range at the same time anisotropy is minimized. Nevertheless, n 
or e is always finite without approaching either zero or infinity. It is the result of crushing 
the object to a conducting plane instead of a line so that there is no singular point in 
the coordinate transform. The small anisotropy together with the finite range of n makes 
realization of the cloak easier. 

To test the effectiveness of the designed cloak, suppose the cloak is Afim by 1.5/im, i.e. 
the length unit in the forehead discussion is in fim and the n profile given before is defined 
relative to silica glass {Si02) of e^ef = 2.25. We ignore the anisotropy in the cloak and only 
keep the part of permittivity. In this case, the permittivity of the cloak varies from around 
1.5 to 4.4. This range of permittivity can be obtained effectively by etching or drilling 
subwavelength holes of different sizes along the direction of E-field in a high dielectric, e.g. 
Si. Outside the cloak, it is again the silica glass as background material. Moreover, the inner 
surface of the cloak is coated by a highly reflective metal. To the observer, it is perceived 
as if this is the actual ground plane. This is relevant to the situation of routing light at our 
own will in optical integrated circuits. Suppose a Gaussian beam at a wavelength of 750nm 
(i.e. 500 nm in Si02) is launched at an angle of 45 degrees to the object together with the 
cloak, the total E-field pattern (real part) obtained from a FDTD simulation is shown in 
Fig. [3]^a). The cloak is within the rectangle in dashed line in the figure. The field outside 
the cloak resembles the field as if we only have a fiat ground plane. A reflected beam at 45 
degrees is clearly seen. Moreover, the field inside the cloak shows the interference pattern 
between the incident and the reflected beam. The pattern is simply squeezed upwards if it 
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FIG. 3: (color online), (a) E-field pattern with the cloak located within the rectangle in dashed 
line when a Gaussian beam is launched at 45 degrees towards the ground plane from the left, (b) 
E-field pattern when only the object (a reflecting surface of the shape of the inner cloak boundary) 
is present without the cloak. All lengths are in /im. The width of the beam is around 4/irn- at 
a wavelength of 750nm. The medium above the ground plane and outside the cloak is Si02 of 
n = 1.5 

is compared to the field when only a flat ground plane is present. The corresponding E-fleld 
pattern with only the object present is also shown in Fig. [3]^b) for comparison. Without 
the cloak, the incident beam is deflected and split into two different angles. Therefore, our 
designed cloak successfully mimics a flat ground plane. 

We emphasize that our strategy is valid in both geometrical and wave optics. In fact, it is 
broadband in nature. Our design works for a wide range of frequencies since the frequency 
dispersion (and also the losses) of the dielectrics can be made very small. For example, we 
can launch a Gaussian wave packet centered at the aforementioned frequency to the object 
and the cloak. The transverse width of the packet is set to 2/im with a duration of 20/s 
(8 periods of central frequency). It results an incident packet with circular spatial proflle at 
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FIG. 4: (color online), (a) E-field pattern of the incident wave packet at time t = 0. A Gaussian 
wave packet with 2fim in spatial width and 20/s in time width (at a central wavelength of 750nm) 
is launched at 45 degrees towards the ground plane from the left, (b) E-field pattern showing the 
reflected packet from the cloak at time t = 37/s. (c) E-field pattern at time t = 37 fs when the 
cloak is absent. All length scales are in fim. The medium above the ground plane and outside the 
cloak is Si02 of n = 1.5. 

time t = 0, as shown in Fig. Il](a). It is traveling at 45 degrees towards the object. After a 
duration of 14.8 periods at t = 37/s, the packet is reflected back at 45 degrees and it remains 
the same size of the Gaussian shape with only a small distortion (or dispersion) as shown 
in Fig. m^b). This is certainly a time-domain cloaking effect. For comparison, we have also 
shown the field pattern at t = 37/s if the cloak is absent in Fig, Hl^c). The incident wave 
packet is split into two elongated packets. There is nearly no distortion for the reflected 
wave packet in Fig. Il](b). However, if a thinner cloak or bigger object is used, it is expected 
that the distortion will become larger. It is due to a larger anisotropy neglected in our cloak 
profile. 
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In conclusion, a cloak is designed to mimic a flat ground plane. The involved coordinate 
transform induces no singular values in the material proflle. Moreover, the quasiconformal 
map is the optimal coordinate transform which can minimize the anisotropy of the physical 
medium. The small anisotropy is further neglected so that the cloak can be synthesized by 
only isotropic dielectrics. It can relieve the loss issue by avoiding the usage of resonating 
elements and is a step closer to do broadband transformation optics in optical frequencies. 
The quasiconformal map should also be a convenient technique for both polarizations as 
well and in other applications in transformation optics. 

This work is supported by the Croucher Foundation fellowship from Hong Kong. 
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